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Abstract 

We determine the expected curvature polynomial of random real projec- 
tive varieties given as the zero set of independent random polynomials with 
Gaussian distribution, whose distribution is invariant under the action of the 
orthogonal group. In particular, the expected Euler characteristic of such ran- 
dom real projective varieties is found. This considerably extends previously 
known results on the number of roots, the volume, and the Euler characteristic 
of the real solution set of random polynomial equations. 
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1 Introduction 

1.1 Real zeros of random polynomials 

The study of the real zeros of random polyomials started with a paper by Bloch and 
Polya [6] who investigated polynomials in one variable with independent random 
coefficients from the set {—1,0, 1}. The question was further investigated by Little- 
wood and Offord [24, 25] who estimated the average number of roots with respect to 
different probability distributions of the coefficients. The first asymptotically sharp 
estimate on the average number of real roots was obtained by Kac [19]. He showed 
that the expected number of real roots of a random degree d polynomial is asymp- 
totically f In d if the coefficients are independent and standard normal distributed. 
Erdos and Offord [12] obtained similar results for random coefficients in { — 1,1}. 
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Maslova [27, 28] proved that Kac's asymptotic result in fact holds for large classes 
of distributions. For more details we refer to the textbook by Bharucha-Reid and 
Sambandham [5]. Edelman and Kostlan [11] give a very nice account of this and 
related work. 

The paper by Shub and Smale [35] was a breakthrough for the study of real 
roots of random systems of polynomial equations. The key point in that work is 
an assumption on the underlying probability measure that is very natural from a 
geometric point of view, namely invariance under the action of the orthogonal group. 
This probability measure was first suggested by Kostlan [22]. By sharp contrast 
with Kac's result [19], the expected number of real roots of a random degree d 
polynomial / = Yla=o fa^'^ turns out to be exactly ^/d if the coefficients /„ are 
independent centered Gaussian random variables with variance (^) . To understand 
the invariance, the polynomials should be interpreted as bivariate forms (having 
roots in P^). The resulting probability distribution on bivariate forms is invariant 
under the action of the orthogonal group 0(2). 

More generally, let Hd^n denote the vector space of homogeneous real polynomials 
of degree d in the variables Xq, . . . , X„. Let / = faXg" ■ ■ ■ G Hd,n be 
such that the coefficients /„ are independent centered Gaussian random variables 
with variance (^) := ^^|^'^ , . The induced probability distribution on the space 



of forms of degree d can be shown to be invariant under the natural action of the 
orthogonal group 0(n + 1). We will say that such / is a Kostlan distributed random 
polynomial. Gonsider a system fi{x) = 0, ...,/„ (x) = of Kostlan distributed 
random polynomials. Shub and Smale [35] showed that its expected number of 
real roots equals ^/dl■ ■ ■ dn, i.e., the square root of the product of the degrees of 
the polynomials. This result was also found by Kostlan [22] in the case where all 
polynomials have the same degree. In fact, Shub and Smale's result was a byproduct 
of a probabilistic analysis of nonlinear condition numbers that control the cost of a 
projective homotopy method to solve systems of polynomial equations, see also [7]. 
We remark that the above choice of invariant probability measure seems also natural 
from the point of view of physics [8] . 

The results by Kostlan [22] and Shub and Smale [35] on the expected number 
of real roots in the setting of an invariant probability measure have been extended 
to multihomogeneous systems by McLennan [29] and, partially, to sparse systems 
by Rojas [33] and Malajovich and Rojas [26]. The work of Kostlan [23] contains a 
classification of Gaussian invariant random polynomials along with further results. 
Recently, Azai's and Wschebor [3] gave a new proof of the Shub-Smale theorem based 
on the Rice formula from the theory of random fields. Wschebor [43], for the first 
time, analyzed the variance of the number of real roots. 

1.2 Underdetermined random polynomial systems 

Gonsiderably less is known for the underdetermined case fi{x) = 0, . . . , fsix) = 
(s < n) where the set of solutions is a real algebraic variety of positive dimension. 
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As a measure of its size different choices come to mind. One possible choice is 
the volume, which is finite when the solution set is interpreted in projective space. 
Another generalization of cardinality to higher dimensional solutions sets is the 
Euler characteristic. (This generalization is not only natural from the topological, 
but also from the computational complexity point of view, as shown in [9].) Both 
of these measures have been considered already. 

Kostlan [22] showed that for a system of Kostlan distributed /j, the volume of 
the projective solution set has the expectation 

E vol(Z(/i, ...,/,)) = vol(P"-^) (1) 

in the case di = . . . = ds = d. (Unfortunately, Kostlan never published the proof. 
A proof for the more general case with possibly different d; has been given in [32].) 

Podkorytov [31] considered any centered Gaussian random polynomial / that 
is invariant under the action of the orthogonal group and determined the expected 
Euler characteristic of its zero set 2{f) in projective space P'*. Podkorytov showed 

Ex{2{f)) = In{VS)/In{l) for odd n, (2) 

where 6 is the parameter of / (see Definition 4.5) and In{^/S) := /o'^(l — x'^)^2~dx. 
If / is Kostlan distributed, then its parameter equals the degree d. If n is even, 
2[f) is almost surely a compact odd-dimensional manifold and therefore its Euler 
characteristic vanishes. Apparently, Podkorytov was not aware of Shub and Smale's 
work [35]. His proof is based on some tricky application of Morse theory. An 
application of Podkorytov's result can be found in [20] . 

1.3 Main results 

The original motivation of the present work was to extend Podkorytov's result [31] 
for hypersurfaces to projective varieties of higher codimension. However, a direct 
use of Morse theory as in [31] does not seem feasible. It turned out to be essential 
to study the more general notions of curvature coefficients and curvature polyno- 
mial of a projective variety. Our main result (Theorem 1.1) considerably extends 
and unifies the previously known results on the number of roots, volume, and Euler 
characteristic. It determines the expectation of the curvature polynomial of a ran- 
dom projective variety under an invariant probability measure. Before stating our 
result, we need to explain the notion of curvature coefficients. 

In a seminal work, Weyl [41] derived a formula for the volume of the tube 
r(M, a) := {y G 5" | dist(y, M) < a} of radius a around an m-dimensional compact 
smooth submanifold M of the sphere 5". Let s := n — m denote the codimension 
of M. Weyl [41] proved that, for sufficiently small a > 0, vol(r(M, a)) can be 
written as a linear combination 

vol(r(M,a))= J2 Ks+e{M)Jn,s+e{a) (3) 
0<e<m, e even 
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of the functions Jn,k{'^) '■= Jq (sin p)'^~^ (cos p)'^~'' dp. The coefficients Ks-\-e{M) 
depend on the curvature of M in and will be thus called the curvature coefficients 
of the submanifold M. The functions Jn,k{oi) determine the volume of tubes around 
namely vol(^(5"-^ a)) = J„,fc(a). 
Let On-i ■= 27r"/^/r(n/2) denote the (n — 1) -dimensional volume of the unit 
sphere S^"^. Following Nijenhuis [30] we rescale the curvature coefficients by 

Me(M) := ^ Ks+e{M) for < e < m, e even (4) 

and define the curvature polynomial of the compact smooth submanifold M of 5" 
by setting 

fi{M;T):= Yl /^e(M)^^ (5) 

0<e<m, e even 

For example, a subsphere 5™ of 5" satisfies p{S'^]T) = 1. One can show that 
the constant term of the curvature polynomial describes the volume of M, namely 
MM) = 0-'yo\{M). 

The generalized Gauss-Bonnet theorem of AUendoerfer and Weil [2] and Her- 
glotz [16] implies that the Euler characteristic xi^) of ^ can be retrieved by 
evaluating the curvature polynomial of M at 1: we have x{M) = 2 p{M; 1) if m is 
even, cf. Theorem 2.1. 

Using the canonical isometric 2-covering map n: ^ P", we define the curva- 
ture polynomial of a compact smooth submanifold M of real projective space P" by 
l^{M;T) := pii:'^ {M);T)) . It is then easy to see that Pq{M) = vol(P™)- Wol(M) 
and x(M) = n{M- 1). 

Suppose now that / = (/i, . . . , fs) G H^^^n x • • • x Hd^,n is a Gaussian random 
system of polynomials. It can be deduced from Sard's lemma that the hypersurfaces 
Z{fi) intersect transversally almost surely, in which case the real projective zero set 
. . . , /s) C P" is a smooth projective variety of pure codimension s, or empty. 
To avoid this case distinction we define /i(0;T) := 0. 

We can now state our main result. 

Theorem 1.1 Suppose that fi € Hii^^ni ■ ■ ■ i fs G ^ds,n < n) arc independent 
centered Gaussian random polynomials with a distribution that is invariant under 
the action of the orthogonal group 0{n + 1). Let 6^ denote the parameter of fa- 
Then the expected curvature polynomial of the projective variety Z{fi, . . . , fg) in P" 
is determined by 

All 

EM2(/i,...,/.);r) = n _ modT-+\ 
In particular, E He{Z{fi, . . . , fs)) depends only on 6i, . . . , 6s and e. When all pa- 
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rameters 5^ are equal to S, we obtain 



A;=0 

where cj,'^ = 1 and C^^^ = sis+2)(s+4ywk~2) ^ ^ ^ 

For instance, the theorem imphes Evol(^(/i, . . . , fs)) = ((^i • • • ^sY^"^ vo^P"^'*), 
which generaUzes (1). In the case where aU parameters equal 5 and ra — s is even, 
the expected Euler characteristic satisfies 

n — s 

Ex(^(/i,...,/.)) = 5^'^Y.^\:\l-5f 

k=0 

= (-1)^ eg, 5t+0(,5t-i) (<5^oo). 

2 

1.4 Methods of proof 

The proof of our main theorem is inspired by Azais and Wschebor's [3] new proof of 
the Shub-Smale theorem based on the Rice formula from the theory of random fields. 
Starting from Weyl's tube formula (3), we derive a version of a "Rice formula" for 
curvature coefficients (Theorem 5.2) and proceed by a probabilistic analysis of that 
formula, making heavily use of the orthogonal invariance. 

In fact, for the proof of Theorem 1.1, it is sufficient to consider the case of one 
equation. This follows by employing a version of Chern's [10] kinematic formula 
of integral geometry for real projective space P". According to Nijenhuis [30], the 
expected curvature polynomial of the intersection of a compact smooth submani- 
fold M of P" with a randomly moving compact submanifold of P" is a truncated 
product of the curvature polynomials of M and N, cf. Theorem 2.2. Actually, this 
was shown by Chern and Nijenhuis for submanifolds of Euclidean space, but it also 
holds for projective space, cf. Santalo [34]. Details on this can be found in the 
monograph by Howard [18]. 

Despite the considerable simplification of the arguments for the case of a hyper- 
surface, we develop the proof via Weyl's tube formula and the Rice formula in full 
generality. Besides some intrinsic interest in obtaining a self-contained probabilis- 
tic proof, the main reason for doing so is that this avenue will also allow to treat 
higher moments (variance), as recently done so by Wschebor [43] for the case of a 
zero dimensional solution set. We thus lay the ground for a planned future paper 
which will investigate under which conditions the curvature polynomial (or Euler 
characteristic) is well approximated by its expectation. 

We also remark that Theorem 1.1 can be quickly derived from the knowledge 
of the expected Euler characteristic of a random projective hypersurface 2{f), for 
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an invariant Gaussian /, as derived by Podkorytov [31], cf. (2). This reduction — 
which wc found first — is again based on the kinematic formula and the generahzed 
Gauss-Bonnet theorem. We present it in §6. 

Finally, we remark that there is some connection of our work to the the study of 
the geometric properties of random fields. Indeed, a random polyomial /(Xq, . . . , X^) 
can be seen as a "polynomial random field" defined on or on the sphere S"". 

A central topic in Adler's book [1] is the study of the Euler characteristic (and its 
variation called IG characteristic) of the "excursion sets" {x & D \ f(x) >u} where 
/ is a real valued Gaussian random field on R"'+-'^ , it G M, and D C R""*"-^ a compact 
domain with smooth boundary. Part of the interest comes from the insight that the 
expected Euler characteristic is a useful approximation to the distribution of the 
maximum of /, cf. [1]. Worsley [42] describes some applications of the statistics of 
the Euler characteristic of excursion sets in to astrophysics and medicine. Adler 
and Taylor [40] have considerably extended and unified the previous results on the 
expected Euler characteristic of excursion sets to the general framework of a cen- 
tered regular Gausssian field on a compact manifold. We realized that it is possible 
to deduce the expected Euler characteristic of a random projective hypersurface, 
that is, Podkorytov's result [31], from the general result in [40, Theorem 4.1], cf. 
Remark 6.4. We remark that our proof is methodically quite different from the ones 
by Podkorytov, and Adler and Taylor. Both of them rely on Morse theory, while we 
analyze expected tube volumes. As already mentioned before, the direct application 
of Morse theory does not seem feasible for the probabilistic analysis of projective 
varieties of higher co dimension. 

The structure of the paper is roughly as follows: In §2 we recall relevant facts 
from differential and integral geometry. Sections 3-4 prepare for the proof of The- 
orem 1.1, which is then given in §5. Hereby, §3 develops the necessary facts about 
Gaussian random vectors and symmetric matrices that arc invariant under the action 
of the orthogonal group. In §4 we give a discussion on invariant random polynomi- 
als. In §6 we show how to quickly derive Theorem 1.1 from the knowledge of the 
expected Euler characteristic of a random projective hypersurface for an invariant 
centered Gaussian random polynomial. 
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2 Background from differential and integral geometry 



2.1 Weyl's tube formula and curvature polynomials in spheres 

For the following material from differential geometry we refer e.g. to [21] or [36, 37]. 
Let M be a compact smooth m-dimensional submanifold of 5", interpreted as a 
Riemannian submanifold. We denote by T^M the tangent space of M at a point 
X G Af and write (TxM)-^ for its orthogonal complement in TxS"". The curvature 
of M at X is described by the second fundamental form of M at x, which is a 
trilinear map 11m {x): T^M x T^M x {T^M)^ M that is symmetric in the first 
two components. In terms of local coordinates, IIm can be described as follows: let 
u = {ui, . . . ,Um) I—*- ipiui, ■ ■ ■ ,Um) G M C be a local parametrization of M. 

Then, for any unit normal vector ly G (T^M)-^, 

llM{x){da,d,3,u) = {dli3(p,u), (6) 

using the shorthand notation da ■= and d'^ ^ := d^^^^^. By the Weingarten map 
of M at x in direction v we understand the self adjoint linear map Lm{x, v) ■ T^M — ^ 
TxM characterized by 

(Lm(x, u){V), W) = 11m{x){V, W, v) for V,W ^ T^M. (7) 

In a seminal work, Weyl [41] determined the volume of the tube 

T{M, a) := {y G 5" | dist(i/, M) < a} 

around M for a sufficiently small radius a > 0. He proved that 

v„i(T(M,„))=/ ff ''-'•';'f-l'^-^^'''^^ is.(.),,mi.). (8) 

where a = tana, Sx denotes the unit sphere in (TxM)-^, and s = n — m is the 
codimcnsion of AI in S*" (see also [15]). Moreoever, Weyl showed that the tube 
volume vol(r(M, a)) can be written as a linear combination of the linearly inde- 
pendent functions Jn,s+e{oi) with real coefficients Ks+e{M), for < e < m, e even, 
cf. Equation (3). The lowest order coefficient Ks{M) equals Os_ivol(M), which is 
intuitively plausible. 

We will call i^s_|_e(-M) the curvature coefficients of the submanifold M of S"". In 
order to justify the naming of these coefficients, we remark that, after some rescaling, 
ifs_l_e(M) is an isometric invariant of the Riemannian submanifold M oi S'^. More 
precisely, for e > 0, 

s(s + 2)---(s + e-2) (• 

t-^s-l Jm 

with some function k(,: M ^ M whose value at x £ M depends only on the difference 
of the Riemann tensor of M and the Riemann tensor of 5" restricted to M, at x, 
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cf. [41]. (We will not need this observation in the following.) It is important to realize 
that these curvature coefficients are not "absolute" invariants of the Riemannian 
manifold M. One can show that for a subsphere M = S"* of we have Ks+eiS"^) = 
for e / 0. 

For our purposes, it will be more useful to rescale the curvature coefficients 
as done by Nijenhuis [30]: we define the curvature polynomial fi{AI;T) as in the 
introduction (Equations (4) and (5)). It follows from the above that the constant 
term of the curvature polynomial describes the volume of M: we have Ho{M) = 
O-^ vol(M). We note that niS"^; T) = 1 for a subsphere S"" of 5". 

The curvature polynomial also encodes the Euler characteristic of M in a simple 
way. The following statement can be deduced from the generalized Gauss-Bonnet 
theorem of Allendoerfer and Weil [2] and Herglotz [16]. We provide the proof in the 
appendix. 

Theorem 2.1 Let M be a compact smooth suhmanifold of 5" of even dimension m. 
Then the Euler characteristic of M can be expressed as x{^) — 2 l^iM; 1). 

2.2 Principal kinematic formula of integral geometry for spheres 

One of the main goals of integral geometry is to compute integrals of the form 
Jq I{Mr\gN)dg, where M and A'' are compact smooth submanifolds of a homogenous 
space with respect to the action of a Lie group G, I is some integral invariant, and 
the integration is with respect to the invariant measure on G. Kinematic formulas 
provide answers to this question in the form I{M fl gN)dg = CkIk{M)Jk{N) 
with integral invariants I/., related to /. For a comprehensive treatment of this 
subject we refer to Santalo's book [34]. A unified treatment of kinematic formulas 
in homogeneous spaces has been given by Howard [18]. 

Chern [10] and Federer [13] proved a general kinematic formula for submanifolds 
of Euclidean space with respect to the group of motions. Nijenhuis [30] pointed 
out a particular elegant formulation of this kinematic formula. He observed that, 
after some rescaling of integral invariants, CkIk{M)Jk{N) can be interpreted as a 
reduced polynomial multiplication. This leads to a great deal of simplification in our 
calculations, as the formulas for the coefficients Cfc turn out to be quite complicated. 

For submanifolds of the sphere and the orthogonal group, the kinematic formula 
takes exactly the same form as for Euclidean space. An indication of this at first 
glance astonishing fact can be found, somewhat hidden, in Santalo [34, IV. 18.3. 
p. 320] for the special case of the intersection of domains. Howard [18] clarified 
this phenomenon by establishing a general transfer theorem according to which the 
Chern-Federer kinematic formulas hold in all simply connected homogeneous spaces 
of constant sectional curvature and not just in Euclidean space. 

The kinematic formula for spheres allows the following beautiful formulation. 
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Theorem 2.2 Let M and N be compact smooth submanifolds of S'^ of the dimen- 
sions m and p, respectively, such that m+p>n. Then we have 

J ii{M n gN- T)dg = ii{M- T)n{N; T) mod 

where the integration is with respect to the Haar measure on the orthogonal group 
0{n + 1) scaled such that the volume of 0{n + 1) equals 1. In particular, we have 



I- 



lx{M n gSP; T)dg = n{M; T) mod t"^+p-"+i. 
We note that Poincare's formula 

0-Vp_„vol(M n gN)dg = O-Vol(M) 0;Vol(iV), 



is a special case of Theorem 2.2, obtained by comparing constant coefficients. 

Weyl's tube formula and the kinematic formula immediately extend from <S" to 
the real projective space P", using the canonical isometric 2-covering map tt: 5" — 
P". We define the curvature polynomial of a compact submanifold M of real pro- 
jective space P" by 

fi{M;T) := fi{7T-\M);T). (9) 
This definition gives the appropriate scaling, as 

Mo(M) = M^-\M)) = ±-yol{Tr-\M)) = ^vol(M) = _^vol(M). 

Moreover, it is easy to see that the kinematic formula of Theorem 2.2 also holds 
for P". Finally, as tt : 5" — P" is a two sheeted covering map, Theorem 2.1 implies 
that x(M) = ^x(7r-i(M)) = fiin-^M); 1) = ^(M; 1). 

Remark 2.3 It is possible to derive Theorem 2.2 from the kinematic formulas given 
in Santalo [34, (15.72), p. 269 and IV.18.3. p. 320] for the intersection of domains 
in Euclidean space and spheres, respectively. Santalo [34, p. 222, p. 302] assigns 
to a compact domain Q in or 5", bounded by a smooth hypersurface dQ, the 
following integral of mean curvature 

Mf«"(5Q):= (^~^^ l^^ai{Ki,...,Kn-i)d{dQ). 

Hereby, (Ti{ni, . . . ,Hn-i) stands for the ith elementary symmetric function in the 
principal curvatures Kj of the hypersurface dQ. The curvature coefficients ifi+i(M) 
of a a compact submanifold M of of codimension s can be related to Santalo's 
integral of mean curvature of the tubes around M as follows 

i^,+i(M) = f"" T ^"i hm Mf-(ar(M, a)). (10) 

\ I I a— >0 
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Herebye, Kj{M) is defined by (3) if j > s and j — s is even. Otherwise, we set 
Kj{M) = 0. The proof is as in [39, V§4], where a similar result is shown for 
Euclidean space. 

3 Invariant Gaussian vectors and matrices 

Here we develop facts about invariant random matrices that will be used in §5 for 
the proof of the main result. 

A random vector is called centered iff its expectation is zero. We call two ran- 
dom vectors equivalent if they have the same distribution. In particular, equivalent 
random vectors have the same expectation and covariancc matrix. We shall write 
X ~ A^(0, cr^) to indicate that X is a real valued Gaussian variable with mean zero 
and variance . 

3.1 Invariant random vectors 

A random vector F G is called 0{n) -invariant iff gV is equivalent to V for all 
5 G 0{n). For simplicity, we assume that V has a density. 

Lemma 3.1 Let V G M" he an 0(n) -invariant random vector. Then V/\\V\\ is 
uniformly distributed in S"'~^ and independent of \\V\\. 

Proof. By invariance, the density of V depends only on \\V\\. □ 

Invariant Gaussian random vectors are easy to characterize by the following well 
known fact. 

Lemma 3.2 A Gaussian random vector V = (Vi, . . . , Vn) is 0{n) -invariant iff it is 
centered and Vi, . . . ,Vn are independent and have the same variance. 

Proof. Suppose V is 0(n)-invariant. Then —V is equivalent to V (take g = —In), 
hence KV = 0. By assumption, gVV-'^g^ is equivalent to VV^ , hence the covariance 
matrix A := K(VV'^) satisfies gAg^ = A for all g G 0{n). A straightforward 
calculation shows that A must be a multiple of the unit matrix (take for g 
rotations in two dimensional coordinate subspaces). The converse follows from the 
fact that Gaussian random variables are characterized by their expectation and 
covariance matrix. □ 

3.2 Invariant random symmetric matrices 

Let S„ denote the space of real symmetric n by n matrices. The Frobenius norm of 
G Sn is defined as \\W\\f ■= i^ij W^jY^^. We wih assume n > 1. 
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Definition 3.3 The parameter S{W) of a random matrix G S„ is defined as 

6{W) := (E {tvWf - E \\W\\l) . 

A random matrix G S„ is called 0{n) -invariant iff gWg^ is equivalent to W for 
all g e 0{n). 

The following proposition classifies all invariant Gaussian W G S„, compare [31, 
Prop. 4.1]. 

Proposition 3.4 Suppose T G S„ with independent Tij and Tij ~ A''(0, 1) for i ^ j 
and Tii ~ N{0, 2), n > 1. Moreover, let r, s G R and Z ~ A^(0, 1) be independent 
ofT. Then W = rZIn + sT is 0{n) -invariant and has the parameter 6{W) = r'^ — s'^. 
Moreover, any 0{n) -invariant Gaussian G E„ is equivalent to one of this form, 
in particular EW = 0. 

Proof. In order to see that T is 0(n)-invariant it is sufficient to check that 
T = gTg^ is equivalent to T for a set of generators g. Since T and T are both cen- 
tered Gaussian it is sufficient to check that they have the same covariancc matrix. 
The group 0{n) is generated by diag(— 1, 1, . . . , 1) and the rotations in two dimen- 
sional coordinate spaces. Invar iance under the latter is verified by a straightforward 
calculation of covariances. Moreover, a calculation shows that 

E(trPF)2 = nV -\- 2ns^, E WWfp = nr'^ n{n + l)s^. 

This implies = — s^. 

Suppose now that W is 0(n) -invariant. Let g be the prodiict of the permutation 
matrix P^^-i and diag(ei, . . . , e„) where Si = ±1. Then gWg^ = (sjej W^^jj 
is equivalent to W. Taking vr = id, ei = —1,62 = ■■■^n = 1 we conclude that 
{Wii, W12) is equivalent to {Wn, -Wu), hence EWu = and E (1^111^12) = 0. By 
choosing appropriate tt and Si one can similarly show that E {WijWki) = except 
in the cases where i = j and k = £ (details are left to the reader). Similarly, one 
shows that W is centered. Moreover, by conjugating with permutation matrices, we 
see that 

EWi = E W^^, E = E =■■ s\ E (WuWjj) = E {WUW22) =: {i + j). 
If we can prove that 

EVFfi = E {WUW22) + 2 E Wf^, (11) 

then W has the same covariance matrix as rZIn + sT and we are done. 

For showing (11) suppose without loss of generality n = 2. Put W = gWg^ 

, /COSCP — SinwA ~ 9 TTT ■ ^TTT 

where g = . • Then Wu = Wn cos ip+W22 sm Lp—2W\2 cos 03 sm ld. 

\smp) cost/? y 

Hence 

ETyfi = E Wl^ (cos^ if + sin^ </?) + 2 (2 E WI2 + E (VFn VF22)) cos^ sin^ p, 
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Using that cos^ + sin^ (p = 1 — 2 cos^ (^sin^ (p, the assertion (11) follows. □ 

Lemma 3.5 Suppose W,Wi, . . . ,Ws € S„ are random symmetric matrices and 
Ai, . . . , As G M. Farther, let u he a real random variable. Then: 

(i) S{W + uln) = S{W) +Eu^ + lE{u trW) . 

(ii) S{XiWi + - ■■ + XsWs) = XjS{Wi) + - ■■ + Xl6{Ws) ifWi, ...,Ws are independent. 

Proof. (i) Fut W = W + Then {tiWf = (trVF)^ + 2nutiW + n'^u'^. 

Moreover, ||W||p = + ^utrW + nu'^. The first assertion follows. 

(ii) Put W := ^<TWa- We have {tiWf = J2a,T ^'^^r {tvW„){trWr). By inde- 
pendence and since arc centered we get E (tiWf =J2a^l^ (trVF^)^. Similarly, 
we obtain EW^^ = E((W^^)h)'. Hence E = A^ The sec- 

ond assertion follows. □ 

Here is a further result stating that stochastic independence follows from invari- 
ance. 



Lemma 3.6 Consider a random (n, V, W) G M x x E„ witii joint Gaussian 
distribution such that (u, gV, gWg'^) is equivalent to {u,V,W), for all g G 0{n). 
(We will call such {u, V, W) 0{n)-invariant.) Then: 

(i) V is independent of u and W. 

(ii) If E {utrW) = 0, then u and W are independent. 

Proof. V and W are centered by invariance and we may assume w.l.o.g. that u 
is centered, (i) Taking g = —In we see that {u,—V,W) is equivalent to {u,V,W). 
Hence E {uVi) = and E {ViWjk) = 0. 

(ii) We argue similarly as in the proof of Proposition 3.4. Using that {u,gWg'^) 
is equivalent to {u, W) for g = diag(ei, . . . , e„) with = ±1 we get E (uWij) = for 
i ^ j. By taking for g a permutation matrix we conclude E (uWu) = E (uWn). By 
assumption, E {uWn) = n~^E {uti W) = 0. Hence u is uncorrelated with all Wij. □ 

The following is a consequence of Lemma 3.6 and Gaussian regression. 

Corollary 3.7 Consider a random {u, W) G Mx E„ with joint Gaussian distribution 
such that {u,gWg'^) is equivalent to {u,W), for all g G 0(n). Let Wc denote the 
random matrix W conditioned on u = Q. Then Wc is 0{n)-invariant Gaussian and 
has the following parameter 
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Proof. Consider the random matrix W := W — Xuin, where A := ^ ^"^'^ . 
Then E(ntrVF) = 0. Moreover, {u,W) is 0(n)-invariant Gaussian. According to 
Lemma 3.6, u and W arc independent. Hence Wc has the same distribution as W. 
Moreover, by Lemma 3.5, 

S{W) = d{W) + X^Eu'' -^E{utTW) = 5{W) - A^En^ = d{W) - ^^^""^V^^ 
as claimed. □ 



3.3 Expected determinant of invariant random matrices 

This section is devoted to the proof of the following result. It is stated in [31], but 
the proof is only vaguely sketched there. 

Proposition 3.8 Suppose W e T^n is 0{n)-invariant and Gaussian, n > 1. Then 
we have 

E det(/n + W)=E{1 + 

where X ~ N{0, 1) is a standard normal random variable. (Note that 5{W) may be 
negative.) 

To prepare for the proof we denote by {k G N) 

the A;-th absolute moment of a standard normal random variable X ~ iV(0, 1). In 
particular, 

72„ = l-3-5---(2m-l) = (2m)!/(2"*m!), 7^0^ = 2(27r)'=/2. (13) 

The following lemma about the higher moments of Gaussian vectors is well 
known, cf. [1, p. 108]. 

Lemma 3.9 Let (Yi, . . . , Yn) be a centered Gaussian vector. Then E{Yi ■ ■ ■ 1^) = 
if n is odd. Otherwise, if n = 2m, 

• • • Y2m) = J2^{Yi,Yi,) ■ ■ ■ E(y,,_,y,,^), 

where the sum is over all 72m = (2m) ! / (2'"to!) different ways of grouping Yi , . . . , Y2m 
into m pairs. 

Clearly, the lemma implies that E det = for any centered Gaussian matrix 
W G M'*^" if n is odd. 
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Corollary 3.10 The random matrix T E'En of Proposition 3.4 satisRes E det T = 
(— l)"*72m if n = 2m and E det T = otherwise. 

Proof. We have E detT = sgn(7r) E(ri ,^(1) • ■ ■T2m,TT{2m))- Lemma 3.9 and 
the fact that the entries of T arc independent imply that only the products of trans- 
positions of the form vr = (nji) • • • (imjm) with {hji,..., im,jm} = {1,2,..., 2m} 
contribute to the sum. The contribution of each such tt is (—1)™ and there are 72m 

such TT. □ 

Proof of Proposition 3.8. According to Proposition 3.4 we may assume W = 
rZIn + sT. Then S{W) = — s^. By expanding the determinant it is easy to see 
that 

det(4 + W) = det((l + rZ)I^ + sT) = ^(1 + rZ)"-l-^l sl-^l detr^,^, 

K 

where the sum is over all subsets of {1, 2, . . . , n} and Tk^k denotes the principal 
submatrix obtained by selecting the rows and columns in K. Using Corollary 3.10 
and taking into account the independence of Z and T we obtain 

[n/2\ 

E det(7„ + T^) = ^ Me (1 + rZr-^^ s^^V 1)^72,- ■ 

i=o ^ 

Suppose Y N{0, 1) is independent of Z. Using i = we have (— l)-^72j = 

E(iy)^-' and hence 

Ln/2J 

detiln + W) = J2 (^Mii + rZ)''-''^E{isYf^ 
= E(l + rZ + isF)" 

\n/2\ 

= E [lY{rZ + isYfK 
On the other hand, using that (2^) 72fc72j-2fc = H) l2j for < A; < j, we get 

E {rZ + ^sYf^ = ^ (^^\ r'>^^,, (-l)^"^ s'^-'^ 72,-2. 
fc=o ^ ^ 

I — n V / 



fc=0 

We conclude that 

Ln/2J 



det(7„ + W)=^i '',jE(Vr2-s2z)2j = E (1 + - s^Z)", 
which finishes the proof. □ 



14 



3.4 Two auxiliciry results on random matrices 

Let A G M*^" be a matrix of rank s < n. We denote by A^n G M" the m-th 
row of A and by the orthogonal projection of Aj^ onto the space spanned by 
Ai, . . . , Am-i- Let vol(yli, . . . , Am) denote the volume of the parallelepiped in 
spanned by Ai, . . .,Am. Clearly, vol(^i, . . .,Am) = vol(^i, . . .,Am-i) H^mll- 

Lemma 3.11 Suppose that A £ M*^" is a random matrix with independent stan- 
dard Gaussian entries. Then we have for 1 < m < s and < j < L^^^^J 

■jg A A \ ' \ ^n—m lfnlfn—m+l—2j 



,(vol(Ai,...,^^_l)-— ^2Pt) = 



\\Air^-'J 2(27r)^ 7n-m+i ■ 

Proof. The invariance under rotations of the standard normal distribution shows 
that, conditioning on Ai, . . . , A^-i, the random variable has a standard normal 
distribution in the orthogonal complement of Ai, . . . , Am-i in M", which is of di- 
mension n — m+1 with probability one. Hence E . . . , = Kn-m+i-, 
where we have written K^i := E ||X|| for a standard normal X in W^. An elementary 
computation gives (cf. [3]) 

Put vol„(^) := vol(^i, . . .,Am). Then voU(^) = voU_i(^) P^ll- Hence 

EvoU(^) =e(voWi(^) E (^\\Ai\\/Ai,...,Am-i)) =EvoU_i(^) Kn-m+l, 
which implies 

EvolmiA) = Kn-m+lKn-m+2 ■ ■ ■ = —■ (14) 

Using polar coordinates, we get for almost all values of Ai, . . . , A^^-i 

^A^(\\Air'^/A„...,A^A = r r"-'^^'e-'4 dr 

1 O 

~ T,~ , n-m 7n-m-2j+l- 
^ (27r) 2 

The claim follows by combining this with (14). □ 

The Moore- Penrose inverse of a matrix G M*^" of rank s is defined by ^1^ := 

{AA^)~^ G IR"^^^, cf. [4]. It is characterized by the following properties: AJ^ = Ig 
and A'^A is the orthogonal projection onto the orthogonal complement (kcr A)-^ of 
the kernel of A. Note that (ker^)-*" is generated by the rows of A. Let us denote by 
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A the restriction of the hnear map A: R" — ^ W to (ker^)-*-. Then the hnear map 
At : R« ^ M" is the inverse of A. We note that 



\detA\ = vol(Ai,...,yl,) = -^det(AA^). (15) 

(Proof: A^Ci = Ai, hence Idet^-^l = vol{Ai, . . . , Ag). It is well known that the 
latter equals Vdet AA^.) 

We denote by (A^i, . . . , Ng) the orthonormal basis of (ker A)-^ obtained from the 
rows of A by Gram- Schmidt orthogonalization. This defines the orthogonal map 

A' 

Qa : ^ (ker A)-^, u^N^. (16) 

By composing Qa with the adjoint (A^)^ of the Moore-Penrose inverse we obtain 
the linear endomorphism {A'^Y^Qa'- K'' W. 

The following is the main result of this subsection. It will be used in the proof 
of Theorem 1.1. 

Proposition 3.12 Suppose that A G M*^" is a random matrix with independent 
standard Gaussian entries. Let u G S^~^ be a uniformly distributed random unit 
vector which is independent of A. Then the random variable 

(VdeiAA^, {A^fQAiu)) 

with values in M.xW has the same distribution as 

(vol(Ai,...,A,_i)- pf||,^^-«;) 

where w is uniformly distributed in the sphere S^~^ and independent of A. 

We give some preparations for the proof. Let A G M*^" be of rank s. Pick 
y G (ker A)-*- and define A^^^ as the restriction of A to the orthogonal complement 
of My ® ker A. This can be described by the following commutative diagram (where 
upgoing arrows denote inclusions): 

(kerA)^ W 
T T ■ 

(MyekerA)^ ^ imA^f) 

Suppose now \\y\\ = 1 and let {Ay)-^ denote the orthogonal projection of Ay onto 
the image of A^y\ Then it is easy to see that 

\detA\ = \detA^y^\-\\{Ay)^\\. (17) 
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Lemma 3.13 Let A G R*^" be of rank s < n and y G (kerA)-*- be a unit vector. 
Then the length of the vector v = {A^)'^y G W can be described as 



V = 



det A^y'^ 



detA 



1 



II (^y) 



±1 



Proof. Wc have v'^A = A = y^, hence A^v = y. Thus v can be interpreted 
as the unique solution of the system of hnear equations 



viAi-\ h VsAs 



y- 



Cramer's rule implies that 



\Vi\ = 



vol(^i, . . . , . ..,As 



vol(Ai,...,^) 

Let Ai denote the projection of Ai onto the orthogonal complement of My. Then 

vol(Ai, . . . , A-i, y, A+i, ■■■,As) = vol(Ai, . . . , ^i+i, ...,As). 
It is now sufficient to prove that (cf. (15)) 

s 

J] vol(]4i, . . . , A+i, . . . ,Asf = (det A(«))2. 

i=l 



(18) 



By an orthogonal transformation, we may assume without loss of generality that 
kerA = Ox M""'' and y = (0, . . . , 0^1, 0, . . . , 0) (with 1 at position s). Write A = 
{ttij). Then aij = for j > s and A is given by the matrix iaij)ij<s with respect 
to the canonical bases. Moreover, A^^) is given by the matrix M : — (oij)i<sj<s iii 
]^sx(s-i) rpj^g relation of Binet-Cauchy [4] states that 

s 

^{det Mif = det(M^M), 

i=l 

where Mj stands for the square matrix obtained from M by deleting the i-th row. 
This is exactly the asserted Equation (18). □ 



Proof of Proposition 3.12. Consider the following random process. First 
choose A G M*^" at random with independent standard Gaussian entries. Then 
pick a unit vector y G (kerA)-*" uniformly at random. Then the resulting random 
variable {V det AA'^,{A^)'^{y)) is equivalent to (Vdet AA^, {A^yQA{u)) . 
Equations (15) and (17) give that 

^Jdet{AAT) = I det I| = I det A^^^ | • || (Ay)^ \\. 
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Lemma 3.13 implies that 



with a random unit vector w & ^. It is clear that {A^)^{y) is 0(s)-invariant. Ac- 
cording to Lemma 3.1, this implies that w is uniformly distributed and independent 
of II (^2/)^ II . 

As in the proof of Lemma 3.11, we sec that ||^^|| conditioned on Ai, . . . ,As-i is 
equivalent to ||-^||, where X is standard normal in M"~*+^. It is not hard to see that 
||(Ay)-'-|| has the same distribution as ||X||. Similarly, one shows that |detA(^^| is 
equivalent to vol(^i, . . . , Ag_i). (Note that if g'l, . . . , gs-i,y is an orthonormal basis 
of (ker A)-^, then | det A^^^ \ = vol(^5i , . . . , Ags-i).) Finally, || {Ay)^ \\ is independent 
of Idet^^?')!. □ 



4 Invariant random polynomials 
4.1 Classification 

We briefly describe the classification of invariant Gaussian polynomials and refer to 
Kostlan [22, 23] for more details. This classification is just for illustration and will 
not be needed for the proof of Theorem 1.1. 

Recall that H^^n denotes the vector space of homogeneous real polynomials of 
degree d in the variables Xq, . . . , The orthogonal group 0(n + 1) operates on 
Hd^n (from the left) in the natural way: for / G H^^n and 5 G 0(n + 1) we set 
{gf){X) := f{g-^X), where X = (Xq, . . . 

Definition 4.1 A random polynomial / G n is called 0(n + 1) -invariant (in- 
variant for short) iff gf is equivalent to / for all all g G 0{n + 1). 

Remark 4.2 Any invariant random polynomial / G ^ is centered if d is odd. 
Moreover, if d is even, it is easy to reduce to the case where / is centered, cf. [23, 
§5.3]. For convenience we will therefore additionally require that / is centered. 

Example 4.3 The most natural example of an invariant random polynomial in 
Hd^n is obtained as follows. Write / G H^^n in the form / = J2a fce^°'^ where the 
sum is over all multiindiccs a G N""*"^ such that |q!| := = d. Assume that the 

coefficients fa are independent with centered Gaussian distribution and variance 
Var(/a) = (^) := ^^r^l- The covariance function x ^,{x,y) ^ 

E(/(x)/(y)) satisfies"' 

E(/(x)/(y)) = ^E{Up)x'^y^ = Q^"?/" = (E^.y,)' = (^,2/)'- 
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and it is thus invariant under the action of 0(n + 1). It follows that / is an invariant 
Gaussian random polynomial. We will say that / is Kostlan distributed [22]. Up to 
a scalar, this random polynomial can be characterized by requiring that / G -f^d,n 
is invariant and centered and has stochastically independent coefficients fa (cf. [22, 
Theorem 4.5]). It is interesting that the normal distribution is enforced by the above 
requirements only. 

It is possible to characterize all invariant centered Gaussian random polynomials 

of Hd^n- For doing so, it is helpful to start with some general observations. 

A centered Gaussian distribution on y = is characterized by its covariance 
matrix. In coordinate free language, this corresponds to the choice of an inner 
product on the dual space V* defined by 

F* X F* ^ R, (A, ^x) ^ E„ {X{v)fi{v)) 

Consider now the general situation of a compact Lie group G operating on a real 
vector space V (in our situation G = 0(n + 1)). It is easy to see that a centered 
Gaussian distribution on V is G-invariant iff the corresponding inner product is 

G-invariant, that is, {gu,gv) = {u,v) for all G G and all u,v €V. 
From the representation theory of compact groups it is known that 

1. On an irreducible G- module W there is a (up to a positive scalar) unique 
G-invariant inner product. 

2. Any two nonisomorphic submodules of V are orthogonal with respect to a 
G-invariant inner product of V. 

(The first statement follows from Schur's lemma; for the second see [17, §27, p. 29].) 

Suppose that V splits into a direct sum V = ©^i^Fj of pairwise nonisomorphic 
irreducible submodules Wj. Choose an invariant inner product ( , )i on each Wj. 
Then, according to the above facts, the invariant inner products on V are of the 
form 

parameterized by ci , . . . , Cm > 0. 

We briefly describe the decomposition of the 0{n + l)-module H^^n into irre- 
ducible submodules, cf. [14, §5.2.3]. Write := Y17=o-^i ^^'^ consider the equiv- 

ariant linear map (p: H,i,n Hdn, f ^ r^A/ arising from the Laplace operator A 
on ]R"+^. It is known that the kernel 'Hd,n is irreducible. The elements of 'Hd,n are 
called harmonic polynomials of degree d. The map cp has the following eigenspace 
decomposition 

ld/2\ 

Hd,n = r^nd-2^,n, (19) 
i=0 

which is thus a decomposition of Hd^n into nonisomorphic irreducible submodules 
{Hd-2i,n corresponds to the eigenvalue 2i{n + 2d — 2i — 1)). We conclude that the 
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invariant inner products on H^^n can be parameterized by [(i/2j positive numbers. 
We refrain from explicitly describing these inner products and refer instead to [23] 
for details. To summarize, we see that the invariant centered Gaussian random 
polynomials of -ff^^n can be parameterized by [(i/2j positive numbers. 

4.2 Parameter of invariant random polynomials 

A polynomial / G i^d^n defines a differentiable map /: — R. We denote by 
Df{x): T^S"" M and D'^f{x): x r^5" ^ M the first and second order 

derivative of / at x G S***. They are characterized by 

/ ( 11^ + AFil ) = /(^) + ^Df{x){V) + ^D^f{x){V, V) + o(A^) (20) 

for y G T-rS'" and A G M going to zero. At the point q = (1, . . . , 0) we can identify 
the tangent space TqS^ = x with R". Clearly, Df{q) is determined by the 
gradient vector {dif{q), dnf{q)). 

Lemma 4.4 D'^f{q) is given by the matrix {dlj{q))i<k,^<n - df{q) In- 

Proof. Consider g{t) := f{il+t^)-^/^il,tV)) = {l+t^)''^/'^ f{l,tV) for F G M" of 
length 1. By definition of the second order derivative we have D'^f{q){V, V) = g"{0), 
cf. (20). A calculation yields /(O) = J2ke dltfi^WkVi - df{q). □ 

The following definition is from [31]. 

Definition 4.5 The parameter 6{f) of an invariant random polynomial / G is 
defined as 

\\xrE\\Df{x)r 

' nE/(.x)2 

where x G R'^^^ \ {0} (this is independent of x by invariance and homogeneity). 

Remark 4.6 1. If / G H^^n is invariant then Df{q) is Gaussian with covariance 
matrix (5(/) E f{qf /„. Thus S{f) = ^^^^f^f^ for any 1 < A; < n. 

2. The parameter of the Kostlan distribution equals the degree d. 

Proof. The first statement follows from Lemma 3.2. For the second just note 
that /(g) = /(<i,o,...,o) and dxj{q) = /(d-i,i,o,...,o)- □ 

Lemma 4.7 Suppose f G n is 0{n + l)-invaTiant and n' < n. Then the restric- 
tion f G H(i^n' of f to R" "'"^ is 0{n' + l)-invariarit and has the same parameter, i.e., 
5{f') = S{f). Moreover, if f is Kostlan distributed, then so is /'. 
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Proof. It is clear that the distribution of /' is invariant with respect to the action 
of 0(n' + 1). In order to see that the parameter remains the same just note that 
restricting / to M" means substituting the variables X^'+i, ■ ■ ■ , by 0. The 
assertion about the Kostlan distribution is obvious. □ 

We show now that the first and second order derivatives of / inherit the invari- 
ance property from /. 

Lemma 4.8 Let f € Hd^n be an invariant Gaussian random polynomial. Then 
{f{q),Df(q),D'^f{q)) G M x M" x S„ is 0{n) -invariant in the sense of Lemma 3.6. 
In particular, Df{q) is independent of f{q) and D'^f{q). 

Proof. Consider for fixed g £ 0{n + 1) the transformed polynomial h := gf, that 
is, h{x) = f{g'^x). By assumption, h has the same distribution as /. This implies 
that the random vector {h,dkh,dl^h) is equivalent to {f,dkf,dlgf). (This can be 
shown by expanding / and 5 in a fixed basis of Hfj,n with random real coefficients.) 
We conclude that {h{q), Dh{q), D^h{q)) = {f{q),gDf{q),gD'^f{q)g^) is equivalent 
^o{f{q).Df{q),D^f{(l)). □ 

The following proposition from [31] says that the parameter of / G equals 
the parameter of D'^f{q), up to a scaling factor. Since the proof in [31] is incomplete, 
we provide a different proof in the appendix. The assumption of a Gaussian random 
polynomial is only made for simplifying the statement and could be replaced by 
suitable regularity conditions. 

Proposition 4.9 Suppose that f G H^^n is an invariant Gaussian random polyno- 
mial. Then: 

(i) E{f{q)trD''f{q)) = -n6if)Ef{q)\ 

(ii) 5{D^f{q))=5{f)Ef{qf. 

The next corollary will be crucial in the proof of Theorem 1.1. 

Corollary 4.10 Suppose that f G -ff^.n is an invariant Gaussian random polyno- 
mial. Let Wc be the random matrix D^f(q) conditioned on f{q) = 0. Then we have 
6iWc)=dif){l-5{f))Ef{qf. 

Proof. Put u = f{q) and W = D'^f{q). Lemma 4.8 implies that {u,W) is 0{n)- 
invariant. Proposition 4.9 yields 5{W) = 5{f) Ku^ and E(ittrVF) = — n(5(/)Eti^. 
The assertion follows now from Corollary 3.7 □ 



21 



5 Random real projective varieties 



We give here the proof of the main Theorem 1.1. Starting from Weyl's tube for- 
mula (8) we present in §5.1 a version of a "Rice formula" for curvature coefficients. 
We then proceed with a probabilistic analysis of that formula, making heavily use 
of the invariance under the orthogonal group. In order to do so, we need all the 
auxilary material on invariant random vectors, matrices, and polynomials that was 
collected in §3-§4, except §4.1. 

We remark that the kinematic formula of integral geometry (Theorem 2.2) would 
allow to reduce to the considerably simpler case of one equation. However, in view of 
a further development of the theory (higher moments) , we will not use the kinematic 
formula here but instead give a self-contained probabilistic proof. 

5.1 A Rice formula for expected curvature coefficients 

In a first step we are going to derive a somewhat more explicit form of Weyl's tube 
formula (8) for the zero set of homogeneous polynomials in S*"". 

Let fi, ■ ■ ■ , fs £ IR[^0! • • • ) ^n] be homogeneous real polynomials of the degrees 
di, . . . ,ds {I < s < n). They define a differentiable map /: 5" — M*. For a point 
a; G ^'^ we denote by Df{x) : T^S"" W and D^f{x) : T^S"" x T^S" ^ W the first 
and second order derivative of / at x. 

In the following we assume that x G 5", f{x) = 0, and rankD/(x) = s. Then, 
locally at x, the zero set Z = Z{f) is a smooth Riemannian submanifold of 5" of 
dimension n — s. The kernel of Df{x) equals the tangent space T^Z of Z at x. We 
denote the inverse of the restriction of Df{x) to the orthogonal complement (TxZ)-^ 
of TxZ by Df{x)^ : — > [TxZ)^ (Moore-Penrose inverse, compare §3.4). 

The following, certainly well known lemma, expresses the second fundamental 
form llzix) of Z at X (cf. §2.1) in terms of the Hessian Hf{x): T^Z x T^Z W, 
which we define as the restriction of D'^f{x) to TxZ x TxZ. Since we could not find 
an appropriate reference, we have included a proof in the appendix. 

Lemma 5.1 Under the above assumptions we have for V,W e TxZ and u G (TxZ)-^ 

{Lz{x, u){V), W) = llz{x){V, W, v) = -{u, Df{x)^Hf{x){V, W)). 

The derivative Df„{x) : TxS'^ M can be identified with a vector in T^S" via the 
inner product on TxS"". Suppose (A''i, . . . , A^^) is the orthonormal basis of (TxZ)-^ 
obtained from {Dfi{x),...,Dfs{x)) by Gram-Schmidt orthogonalization. We use 
the orthogonal map Qf{x) : W (TxZ)-^, u J2t=i ^cr-^o- to describe unit normal 
vectors in (TxZ)-^ by coordinates. We thus define a Weingarten map 

Lf{x,u) := Lz{x,Qf{x){u)) 
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of Z at X in direction parameterized by u G 5* ^ (compare §2.1). According to 
Lemma 5.1, the Weingarten map is explicitly characterized by 

{Lf{x,u)iV),W) = -{Qf{x){u),Df{x)^Hf{x)iV,W)) 

= -{{Df{x)YQf{x){u),Hf{x){V,W)), (21) 

where (D/(x)t)^: (T^Z)-^ W denotes the adjoint map ofDf{x)^ : W (T^Z)-^. 

We suppose now that rank Df{x) = s for all x G Z = 2{f), i.e., the hypersur- 
faces 2(fi) intersect transversally. By Sard's lemma, this is the case for almost all /. 
Then Z is either empty or a compact smooth submanifold of 5" of dimension n — s 
and we assume the latter. 

It will be convenient to introduce the following function associated with / 

n.-^-xM^^M ..(rf)- det(id - ||t||L/(x, t/||t||)) 
gf.S xR ^R, gf{x,t) .- ^ ||^||2)(„+i)/2 " 

By the transformation theorem, Weyl's formula (8) for the volume of tubes can 
be concisely rewritten as vol(T(Z, a)) = J^^b 9f^i^ ^ -^a)' where a = tana is 
sufficiently small and Ba := {t G R^ \ \\t\\ < a} is the ball of radius a in M*. 
Combining this with (3) we obtain 



gf d{Z xBa)= V Ks+2j{Z) Jn,s+2j{a). (22) 

This expansion is valid for all < a < 7r/2 since the functions on both sides of (22) 
are analytic. 

We can now state the announced Rice formula for expected curvatures, which 
will allow to determine the expectations ¥.fKs+2j{2{f)) of the curvature coefficients. 

Theorem 5.2 Suppose that f = (/i, . . . , /s) G H^^^n x • • • x Ha^^n is a Gaussian 

random vector. Hence f{x) G is a Gaussian random vector for anyx G S*"", and we 
shall denote its density function by pji^^-^ : ^ M. Let u G 5"*^^ he a random unit 
vector uniformly distributed in the sphere and independent of f. DeEne the function 
i/j : [0, oo) X 5" — M by the following conditional expectation for (r, x) G [0, oo) x 



Hr, x) := Pf^,) (0) Ef,u \^^Jdet{Df{x)Df{x)T) det(id - rLf{x, u)) / f{x) = 
By taking a spherical average we define the function 

* : [0, oo) ^ M, ^(r) ■= ^ f V'(r, x) dS"" . 
Then we have for < a < 7r/2 and a = tan a 



j2 Ei^.+2,(^(/)) Jn,s+2j{a) = O^O,.^ £ (1 +7^')t+W2 



j=0 
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Proof. The proof uses similar ideas as in [1, §5.1] and [3, Theorem 1]. 

Fix a > and consider for /j G -ff^. ru 1 ^ ^ ^ the corresponding map 
/: ^ The fibre integral 



Gfiy):= [ gfd{f-\y)xBa 



is well defined for regular values y eW. We thus need to determine (cf. (22)) 

^ EK,+2j{Z{f)) Jn,s+2j{a) = Ef{Gf{0)). (23) 

j=0 

We will apply the coarea formula (or Fubini's theorem for Riemannian mani- 
folds). Recall that the normal Jacobian NJj(a;) := a/ det{Df{x)Df{x)'^) of / at x 
has the following geometric meaning. Suppose x G S*" is a regular point of / and 
consider the restriction (kei Df(x))~^ M.^ of Df{x) to the orthogonal complement 
of kerZ)/(x). Then NJy^(a:) of / equals the absolute value of the determinant of 
this map, cf. (15). The normal Jacobian of the differentiable map of Riemannian 
manifolds 

F:5"x5„^M^(,x,^)^/(x). 

at {x,t) eS"" xBa satisfies N.]F{x,t) = NJ/(x). 

Consider the following integrable function ipd for 5 > 

" ' ' \ otherwise. 
The coarea formula (cf. [18, Appendix] or [39, III. §2]) applied to F yields 

/ Gfiy)dy=[ [ ^s9fdif-\y)xBa)dy 

= [ ^5ff/NJFd(S"xB„). 

JS"XBa 

Dividing by {25Y and taking the expectation over / with respect to the given Gaus- 
sian distribution, we obtain 

/ EfiGf{y))dy = [ E/(<^5 5/NJf) d(S" x B„). (24) 

K^(>) Jye{-S,S)'> Js^xBa [^") 

For fixed (x, t) G S*" x Ba wc can write the integrand Is{x, t) on the right-hand side 
of (24) as an integral over conditional expectations as follows 

= 7^s I Pf{x){y)^f {gf{x,t)mf{x)/ f{x)=y) dy. (25) 
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By continuity, we get 

VimIs{x,t)=pf^,^{0)Ef(^gf{x,t)NJf{x)/f{x) = o). (26) 

The integrand of (25) is a continuous function of (x, t, y) and therefore bounded 
by some constant M on 5"" x Ba x {y i=i W \ \\y\\ < 1}. Hence (25) is as well 
bounded by M for all < < 1. We may therefore apply Lebesgue's Theorem 
and interchange in (24) the integral over {x,t) e x Ba and the limit for 5 — 
obtaining 

(G/(0)) = hm jl- I (G/(y)) dy = lim / h d{S^ x Ba) 

= 1 p/(,)(0)E^(3;(x,i)NJ^(x)//(x) = 0) d(5"xi?J, (27) 

where we have used (26) for the last equality. 

Note that for r > (with uniform random u G S^~^ independent of /) 

^'/(a:)(0)E/,„(5/(a;,™)NJy(x)/ /(x) = O) = ^^^^2)(n+l)/2 ^(^'^)- 

Hence, using polar coordinates i = ru, the right-hand side of (27) can be written as 
Taking into account (23), this completes the proof. □ 



5.2 Expected chciracteristic polynomial of Weingarten map 

In order to prove the main Theorem 1.1 we will evaluate Theorem 5.2 for independent 
Gaussian polynomials fa having invariant distributions. We write da = S{fa) for 
the parameter of We may assume without loss of generality that E/o-(g)^ = 1 
for all 1 < (J < s (scaling does not change the parameter of fa)- Hence faig) 
is standard normal and the joint distribution of f{q) has the density Pf(q){y) = 
(27r)-^/2exp(-i(y2 + . . .^2))^ particular, p;(,)(0) = {2tt)-'/\ 

We proceed by a sequence of intermediate steps. For (r, u) £ [0, oo) x S^^^ and a 
fixed matrix M G M''^'^ of rank s we consider the following conditional expectation 

S{r, u, M) := Ef (det(id - rLf{q, u)) / f{q) = 0, Df{q) = m) . (28) 

Recall the characterization (21) of Lf{q,u) in terms of Df{q) and D'^f{q). Prom 
Lemma 4.8 we know that D^f{q) is independent of Df{q). Hence the above expec- 
tation may be taken with respect to the distribution of D^f{q) conditioned solely 
on the event f{q) = 0. 
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Lemma 5.3 For (r, it) G [0, oo) x 5* ^ and a fixed matrix Df{q) of rank s we set 
v = iv,,..., vsf := diag(5}/^ . . . , 5^2) {Df{q)YQf{q)iu). (29) 
Then we liave (recall (12) 

S{r,u,Dfiq)) = ^ r''l2j(^''~^'^(i2(^-S.)viy. 

Proof. Proposition 3.8 is the key to this result. As usual let TqZ denote the 

kernel of Df{q) and recall that the Hessian H f{q) was defined as the restriction of 
the bilinear map D"^ f(q) to TqZ x TqZ. For the following introduce an orthonormal 
basis adapted to M" = TqZ®{TqT)^ (or observe that §3.2 could have beeen presented 
in a coordinate-free way). 

The matrix of Lf{q,u) is 0{n — s)-invariant and Gaussian. The same is true 
for the random matrix Lf{q, u)condi which is defined as the random matrix Lf{q, u) 
conditioned on f{q) = 0. In order to apply Proposition 3.8 we need to calculate the 
parameter of Lf(q, n)cond- 

We write Hf{q) = {Hfi{q),...,Hfs{q)). By identifying (bi)linear maps with 
their matrices we obtain from Equation (21) that 

s 

Lf{q,u) = -^^HUiq), 
a=l ^"'^ 

where v = (vi, . . . ,Vs) is defined as in (29). Hence Lemma 3.5 implies 

Sl := S{Lf{q,u)cond) = XI ^iHf<r{<l)cond) 



using obvious notation. Lemma 4.7 tells us that the restriction of f„ to TqZ (whose 
distribution is invariant under the orthogonal group of TqZ) has the same pa- 
rameter 5^ as /(J. Corollary 4.10 gives that S{H ffj(q)cond) = <5o-(l ~ S^), hence 
Sl = J2a(^ ~ '^cr)^CT- Proposition 3.8 implies now with X ~ iV(0, 1) 



S{r, u, Df{q)) =K{l+r^5L X^-^ 
Hence, taking into account that 72^ = E X"^^ , we conclude 

E(l + rx/^X)— = f: (^V)"'''^^^^^' 
which proves the lemma. □ 



26 



5.3 Proof of main theorem 

We prove now the following reformulation of Theorem 1.1 for zero sets in spheres. 

Theorem 5.4 Suppose that f„ G Uda^n Siic independent centered random polyno- 
mials with 0{n+ l)-invariant Gaussian distribution of parameter 6^- Consider the 
random zero set 2{f) C 5" where f = {fi, . . . , fg)- Then the expectation of the 
curvature coefRcient Ks+2j{Z{f)) satisfies 

for < j < [^^J , wiiere the coefRcients cj^^ can be characterized by the generating 
function 

(1 - y)-V2 = f;c(^) Y' = 1 + ^Y + ^Y' + ^Y^ + ... . (30) 
fc=o 

We have C^^^ = 1 and = i^^^^^P^ = ^ for k > 0. 

In the case where all S are equal the result simplifies to 

E K,+2j {Z{f)) = 5'/^ (1 - 5y On-s-2jOs+2j-iCf for 0<j< [ VJ . 

(s) 

where the Cj ' are characterized as the coefficients of the power series 



fe=0 



More specifically, we have Cq*^ = 1 and for k > 

^ g(g + 2)(g + 4)---(5 + 2fc-2) 
k }^i2k 

Proof of Theorem 5.4. Put A := diag{y/5^, . . . ,^/5^). Since 5a^^'^Df^{q) 
is standard normal distributed in M" and the are independent, we can write 
Df{q) = AA, where A G R*^" is a random matrix with independent standard 
Gaussian entries. Note that det{Df{q)Df{qf) = 5i---5s dei{AA^) and Df{q)^ = 

Let u G S^~^ be a uniformly distributed random unit vector which is independent 
of A. The function V' introduced in Theorem 5.2 satisfies for r > 

V'(r, q) = (27r)-^/2(<5i • • • J,)V2 (Vdet^^^ S{r, u, AA)) . 



27 



Lemma 5.3 tells us that 

j=o \ J / 

where v € is the image of u under the linear endomorphism 

A{Dfiq)YQfiq){u) = (AYQAiu) 
of W (recall the definition of Qa in (16)). We make the multinomial expansion 

( ^(1 - 5,)v^y = Yi (i) (1 - ^i)" •••(!- • • • ^s'^^- 

Thus we need to compute for G with \v\ = j 

Ea,u (Vd^tAAT vl''' ■ ■ ■ -uf ^) . 

Proposition 3.12 determines the joint distribution of (V det AA^ , v) for random A 
and u. Accordingly, we write with a uniformly distributed w G S^~^ that is inde- 
pendent of A: 



E^,„ (Vdet^^^ vl""' ■ ■ ■ vf^) = Ea,u (vo1(^i, . . . , ^_i) pxpi^i^? 
It is well known that [41] 

"'^'^ V 1 W s(s + 2)...(s + 2j -2) 



By using Lemma 3.11 we obtain 



2*^1 .. . „2!^s 



On-s 



This formula can be considerably simplified. We put C^"^ := for A; G N. 
Claim. We have 

^.r'^ol^o.^.,.. ^- (" - •) (^)e. ( vct . ■ ■ .f) ^ . . . eg.. 

In order to verify this recall first that 

n+l 

'^--r(^)' ^""v^ ^ 2 ^• 
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Prom r{x + 1) = a;r(a;) and r(l/2) = we get 

r(m + 1) = ml, T{m + 1/2) = (m - l/2)(m - 3/2) • • • 1/2^7?. 
Using the above, it is straightforward to check that 

0„7n = 2(27r)«/2, !^ = (27r)^, (2^). ^£z^ = + 2) • • • (. + 2j - 2). 

Moreover, recall from (13) that (2j)! = 2^ j\^2j and note that Cj^^^ = ^f|t- The claim 
follows by simplifying the formula using the above stated equations in a straightfor- 
ward (but tedious) way. 

Combining what we have shown so far we obtain 

aa-i^(r,g) ^ (l-<5,)----(l-5.)-CW...C(VH. 

We apply now Theorem 5.2. From the 0{n + l)-invariancc it follows that 
tp{r,x) = ip{r,q) for all x G S"", hence ^(r) = tp{r,q). Let < a < 7r/2 and 
put a = tana. By substituting r = tanp we obtain 



I 



a ^s+2j-l 

dr = Jn,s+2j{(x). 



lo (l+r2)(n+l)/2 

We conclude with Theorem 5.2 that 

I I 

2^ ?T 'Jn,s+2j[a) 

^ C„_s_2j(-'s+2j-l 

r r^-i*(r) 



-/ 



C»„_,_2,a+2i-l Jo (1 + r2)("+l)/2 



= ((5i---(5,)V2 ^ (l-<5i)'^^---(l-5,)'^»cW---Cg)j„,,H.2H(a). 

IH<L^J 

By comparing the coefficients of the linearly independent functions Jn,k{(^)j the 
stated formula for E i^s-|_2j(^(/)) follows. 

To settle the case where all 6^ arc equal just note (1 — Y)~^/'^ = Yl'k=o ^k^^ 
and (1 - y)-^/2 = Y.t=o C^k^ implies that cf = Zw\=j C-l^ ' ' ' C^l'^ ■ □ 
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6 Alternative proof of the main result 



We show here that Theorem 1.1 can be quickly derived from the knowledge of 
the expected Euler characteristic of a random projective hypersurface 2{f) for an 
invariant centered Gaussian random polynomial /. The key of this reduction is the 
kinematic formula and the generalized Gauss-Bonnet theorem. 

In this section, Z{f) stands for the zero set in P". Suppose / G Hd^2e+i has 
invariant centered Gaussian distribution with parameter 5. We define Xei^) •= 
Ex(Z(/)). Theorem 1.1 in the case of one equation {s = 1) yields 



Xe{S)=S'/'J2ci'\l-6)K (31) 

fe=0 



Taking into account Equation (30), we get a closed form expression for the generating 
function xi^'-iT) of xe{S) ^ follows: 



oo oo oo 

e=o k=o e=k 

(51/2 



(i-r2)(i- (i-(5)r2)i/2- 



(32) 



This formula can also be readily deduced from Podkorytov's result [31], cf. (2). 
Using the kinematic formula we can prove a stability result for E/Xe(.Z(/)). 

Lemma 6.1 Suppose f G H^^n is 0{n + 1) -invariant and n' < n. Then the restric- 
tion f G Hd^n' of f to +^ is 0(n' -I- 1) -invariant and lias the same parameter. For 
0<e<n',e even, we have E fie{Z{f')) =Ene{2{f)). 

Proof. The first part of the statement was already established in Lemma 4.7. By 
Theorem 2.2 (in the version for P") we have for almost all / G Hd,„, as /i(P" ; T) = 1, 

J f,{Z{f) n gF-'-T) dg = /x(^(/); T) mod T"', 

where the integral is with respect to the Haar measure of 0{n + 1) scaled such that 
the volume of 0(n -|- 1) equals 1. Taking the expectation over / and interchanging 
with the integral over g we obtain 

I %(m(^(/) n yP"'; T)) dg = Ef ii{Z{f)-T) mod T^' . 

By the invar iance of the distribution of /, the integrand is independent of g, hence 
the integral equals E/(^(2:(/) n P"'; T)). This expectation equals E//(/i(Z(/'); T)), 
which finishes the proof. □ 



30 



Lemma 6.2 Suppose that f G Hd^n is 0{n + l)-invariant with parameter S. Then 
we have E^o{Z{f)) = xo{f) and for 1 < k < {n — l)/2 

Efi2ki2if)) = Xk{S)-Xk-i{S). 

In particular, E ^2k{^{f)) depends only on k and S and not on the dimension n of 
the ambient space. 

Proof. We denote by /' and /" the restrictions of / to R^'^+^ and M?'', respectively. 
Theorem 2.1 (in the version for P") imphes that 

k k-l 

XkiS) = Y,^M2if')), Xk-i{S) = ;^E/.2^(^(/"))- 

i=0 i=0 
By Lemma 6.1 we have E ^.2i{Z{f)) = E fi2,{Z{f')) = E H2i{Z{f")) for < i < A; - 1 
and ¥a iJL2k{Z{f)) = E/X2fc(-2(/'))- Subtracting the above two equations, we obtain 
Xk{^) ~ Xfc-i(^) = IE ^2k{Z{f)) as claimed. The assertion for = is obvious. □ 

We proceed now with an alternative proof of Theorem 1.1 
Proof of Theorem 1.1. By Lemma 6.2 the formal power series 



M5;r) :=^E/.2fe(2(/))r2 



-<2k 

k=0 



satisfies nid-T) = (1 - T^) x(<5; T). Equation (32) implies that 

E fi{S; T) = ,5^2(1 _ (1 _ <5)r2)-V2. (33) 

Suppose now that G Hda,n are independent random variables with invariant 
centered Gaussian distribution of parameter 6(j for 1 < o" < s < n. Write / := 
(/i, . . . , /s-i). The kinematic formula (Theorem 2.2 in the version for P") tells us 
that for almost all /, fg 



I 



li{Zif)ngZ{f,y,T)dg = fi{Z{fy,T) ix{ZUs);T) mod T 



•n—s+l 



where the integral is over 0{n + 1) with respect to the Haar measure scaled to 1. 
We take the expectation with respect to / and fg- Taking their independence into 
account we get 



/ 



EfJ^^i{Z{f_)^gZ{fs)■T)dg = ¥.^i{Z{f_)■T) ii{Z{fs)-T) mod T 



n—s+l 



By invariance, the integrand does not depend on g and equals E iJ,{Z{fi, . . . , fs);T). 
We conclude by induction that 

EAi(^(/i, . . . ,/.);r) =Etx{Z{h);T) ■ ■ • E/x(.E(/,);r) mod T^'^+K 

Plugging in the explicit expression (33) for E n{Z{fi);T), the desired formula for the 
expectation of the curvature polynomial follows. Finally note that (1 — y)"*/^ = 

Ek=ocl'^y'- □ 
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Remark 6.3 The above reduction to the computation of the expected Euler char- 
acteristic works for any 0(n + l)-invariant distribution of random polynomials. The 
Gaussian assumption is not needed for the reduction. 

Remark 6.4 We briefly outline how (32) can be derived from a general result of 

Taylor and Adler [40]. Suppose that / G H^^n is 0(n+ l)-invariant with parameter S 
and suppose w.l.o.g. that K f{p)'^ = 1 for p € S'". We consider / as a centered unit 
variance Gaussian field on the sphere S'^. Then / defines a Riemannian metric 
on 5" by gp{Xp,Yp) := E{Xpf ■ Ypf) for tangent vectors Xp,Yp e TPS'". It easily 
follows that gp{Xp, Yp) = 5 {Xp, Yp), where ( , ) denotes the scalar product on TpS'^. 
Hence, with respect to this metric, is isometric to the sphere M in M"+^ of radius 
R = VS . Theorem 4.1 of [40] states that 

n 

with "Lipschitz-Killing curvatures" Cj of M and functions pj (u) related to Hermite 
polynomials. Almost surely, S"" H /~^[n, cxd) is a compact domain with smooth 
boundary 5" n f'^{0), which implies xl-S" n /"HO)) = 2x(5'" n f-^[u,oo)) if n is 
odd. The Lipschitz-Killing curvatures of M can be defined via Weyl's formula for 
the volume of the tubes around M in M"+^ 

n 

vol(r(M, r)) = ^ Cj ujn+i-j r"+^-^' 

j=0 

with := 1 and tok = Ok-i/k for k > 0. A straightforward calculation yields 
Cn = ("^^)R^ if 7 = n mod 2 and Cj = otherwise. Hence we obtain for 

odd n 

^ l<j<n,j odd V 3 J 

It is possible to derive Equation (31) from this, but we omit the details. 
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Proof of Theorem 2.1. W.l.o.g. m < n. We suppose first that n is odd. The 
submanifold M is a deformation retract of the tube Tq, := T{M,q) for sufficiently 
small a > 0, hence x{^) = linia^o x(^a)- The generalized Gauss-Bonnet formula 
applied to the domain Ta in S"^ says that (cf. [34, (17.22), p. 303]) 

0<i<n-l,i even '^-l-' « "'9T„ 
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Going to the limit a — and applying Equation (10) we get {s = n — m is odd) 

^X{M) = J2 (puti + l = s + e) 

0<i<n— 1, i even 

E ?^ Ks+e{M)= Yl ^e{M). 

0<e<m, e even 0<e<m, e even 

In the case where n is even, we argue similarly: The generalized Gauss-Bonnet 
formula (cf. [34, (17.21), p. 303] applied to says that 

^0„x(ra) = voKTe,) + V ,^ / ai(Ki,...,K„_i)d(ar„). 

Hence, by taking the limit a — and using (10), we get (note that s is even) 

l<i<n— 1, i odd 

E ?^ 7) ^3+eW= E '"e(M), 

0<e<m, e even 0<e<m, e even 

which finishes the proof. □ 

Proof of Proposition 4.9. The covariance function r(x, y) := E{f{x)f{y)) is 
a polynomial function that is homogeneous of degree d in both sets of variables 
X G R""*"^ and y G R""*"^. Since the distribution of / is invariant under the action 
of the orthogonal group, we have r{gx,gy) = r{x,y) for all g G 0(n + 1) and 
x,y e M""*"^. It follows from invariant theory that r{x,y) is a real polynomial in 
\\x\\^, \\y\\^, {x,y) (e.g., see [38]). Prom the fact that r is bihomogeneous of degree 
{d, d) it is easy to conclude that r has the following form 

Ld/2J 

r(x,y)= J]/5fc||x||2*^||y|p^(x,y)'^-2*^ (A e K). 
fc=o 

We can express the moments of the partial derivatives of / at g by partial 
derivatives of the covariance function r at ((?, g): we have for 1 < i,j,k,i < n 

E {dl^.m m) = dl,.r{q, q), E {d,J{q) dyj{q)) = dl^jiq, q), 

E {dl.,m dlyjiq)) = dt,.,y,yA^, q). (34) 

In order to show the second claim it is convenient to use the abbreviation aij := 
dl.^.f{q), A := {aij)i<i,j<n. By Lemma (4.4) we have W := D'^f{q) = A-df{q)In. 
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We obtain for the parameter of W 

S(W) = , ^ ("e (trWf - E llM^lli) 
n[n — 1) V / 

= — ^ — - fE (tiAf - E Pill - 2(n - l)dE (f(q)tTA) + n(n - l)d^E f(qf) 
n{n — 1) ^ 

= E ^ («-"^-^- - 4) - ^ (/(«) t^^) + f{qf 

On the other hand, 

i i 

In order to prove the second claim, it suffices to check equahty of the above two 
expressions. Since both expressions are hnear in r, it is enough to check this for 

rk{x,y) = ||.t|P'''||i/|P'^(x, y)'^^^'^. 

A tedious but straightforward calculation yields for 1 < i,j < n, i j, 

rk{q,q) = 1, dx,xirk{q,q) = 2A;, dxiyirk{q,q) =d-2k, 
dt,,,y,yT{q, q) = dt,,.y,y.r{q, q) = {d - 2k){d -2k-l). 

Plugging in this in the above expressions we see that indeed 6{W) = 6{f)E f{q)'^. 
The verification of the first claim is similar and a bit simpler. □ 

Proof of Lemma 5.1. By invariancc of the assertion under the orthogonal group 
it is sufficient to verify the claim at the point q := (1,0. .. ,0) and we may also 
assume that TqZ = ker Df{q) = M"~* x 0*. Hence dn-s+rfaiq) = for 1 < cr, r < s. 

Our assumption rankD/(q') = s allows to apply the implicit function theorem. 
There is an open subset U C containing the origin and a differentiable map 

h: R"-^ DU such that /i(0) = and 



^l + \\uf + \\h{uW 

is a local diffeomorphism of U onto an open neighborhood of q in Z. 

We make the following useful convention on indices: a,/3, 7 run in the range 
1, 2, . . . , n — s while a, r, p run in the range 1, 2, . . . , s. A straightforward calculation 
shows that 
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where the last equahty follows from our assumption TqZ = x 0''. A similar 

calculation yields 

By the definition (6) of the second fundamental form we obtain for V, G TqZ = 

R""-" X 0** and u G (T^Z)-^ = 0""* x W 

llz{x){V, W,iy)=J2 ^adlpK{0)VaWfs. (35) 

a,l3,a 

It remains to express pha{0) by partial derivatives of /. 

By taking the derivative of fp{l,u, h{u)) = with respect to Ua we obtain 

dx^fp{l,u,h{u)) +^dx„_,+Jp{l,u,h{u)) ■ daKiu) = 0. 

r 

Differentiating this with respect to up and taking into account that dahr{0) = we 
obtain after a short calculation at v = 

r 

We may write this as 3^^/io-(0) = — '^pma,pd\^ X/jfpi^)-) where {m^^p) denotes the 
matrix of Df{q)K Plugging this into (35) we get 

o- P a,l3 

= {iy,DfiqyHf{q){V,W)), 

which was to be shown. □ 
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